Role of Correlation on Charge Carrier Transport in Organic 

Molecular Semiconductors 



Ajit Kumar Mahapatro and Subhasis Ghosh 
School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India 

Abstract 

We have investigated the charge carrier transport in organic molecular semiconductors. It has 
been found that mobility is a function of electric field and temperature due to hopping conduction. 
Several theoretical models for charge transport in disordered solids have been debated over the role 
of spatial and energetic correlation in these systems and such correlations have been recently shown 
to explain the universal electric field dependence of mobility. We have compared and evaluated 
the applicability of different theoretically proposed models using very simple experimental results 
and based on our extensive analysis, we have found that correlation is important to explain the 
electrical transport in these systems. 
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Recently, considerable debate 



arrier transport in 

organic molecular semiconductors(OMS), which belong to wide class of material, known as 
disordered organic molecular solids(DOMS), because of two reasons: first, their application 
in display devices and second, the fundamental understanding of charge transport in these 
materials. There are several distinguishing features of these solids, (i) they are composed of 
organic molecules, held together loosely by weak van der Waals type intermolecular coupling 
while the intramolecular coupling is strong; (ii) the absence of long range order in these 
disordered materials lead to the localization of the electronic wave function and the formation 
of a broad Gaussian density of states(GDOS) and the most important one is (iii) the carrier 
mobility /i exhibits a nearly universal Poole-Frenkel(PF) behavior^] 

/ i(i 71 ,T) = / i(0,T)exp( 7 (T)v / F) (1) 

where /x(0, T) and j(T) are temperature dependent quantities, known as the zero field mobil- 
ity and the field activation of the mobility, respectively. In DOMS, hopping among the molec- 
ular sites having comparable energies, describes the transport of charge carriers through the 



GDOS of highest occupied molecular orbital ( 
orbital(LUMO). Recently, several theoretical 



iOMO) and/or lowest unoccupied molecular 
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models for charge transport in dis- 



ordered solids have been debated over the role of spatial and energetic correlation in these 
systems and such correlations have been recently shown to explain the universal electric field 
dependence(Eq.l) of mobility. In this paper, we have examined the charge carrier transport 
in OMS and we show for the first time, by means of simple experiments that correlation is 
important in explaining charge carrier transport in OMS. 

For our study we have chosen metal phthalocyanine(MePc) based thin films as the organic 
solids because beside their potential application for plastic based optoelectronic devices and 
transistor, these materials are (i) the most chemically and thermally stable compounds, (ii) 
suitable for organic single layer structures with different metal electrodes without effecting 
the interface properties and finally (iii) reproducibility of experimental data for devices with 
MePc as organic layers, which is a major problem with most other organic solids due to 
their degradation with time and the interaction of organic materials with the different metal 
electrodes, notablly with low work function metal electrodes. MePc based organic single 
layer devices were prepared by sandwiching a thin organic layer of MePc between indium tin 
oxide(ITO) coated glass substrate and thermally evaporated metal electrodes, using vacuum 



deposition technique. We report the experimental investigation on charge carrier transport 
in hole only devices based on met al/MePc /metal structures. Work functions of ITO and Cu 
are 4.75eV and 4.6eV, respectively, very close to the ionization potential (4. 8eV) of MePcQ], 
hence the ITO/MePc and ITO/Cu interfaces behave an Ohmic contact. The temperature 
dependent current- voltage( J- V) characteristics were studied for ITO/MePc/Al structures 
with different MePc layers, choosing ITO as the anode. By properly choosing contacting 
metals(Cu and Al), current injection and transport due to holes have been investigated. Here 
we have chosen copper phthallocyanine(CuPc) and zinc phthalocyanine(ZnPc) as the MePc 
for these investigations. High purity, sublimated CuPc and ZnPc have been procured from 
Aldrich Chemical Co. CuPc(ZnPc) was evaporated on ITO coated glass substrates from a 
resistively heated tungsten boat. Each layer was deposited in a vacuum chamber at a rate of 
5-10A/s. Thickness of the resulting organic films was between 25 to 500nm. The metals for 
the top electrode(Au, Cu and Al) were placed in a resistively heated tungsten filament and 
evaporated in a vacuum chamber. All evaporations were done at a base pressure of 2xl0 -6 
Torr. We have measured the J-V characteristics by continuously varying the bias and also 
by interrupting the bias between each voltage step and in both cases, we have found same 
results. The J-V characteristics for the ITO/MePc/metal structures with different thickness 
of MePc layers were reproducible at all temperatures and showed no hysteresis. 

Fig.l shows the J-V characteristics of ITO/CuPc/Al and ITO/CuPc/Cu at room tem- 
perature. The experiment consisted of two steps. First: the current due to hole injection 
from ITO was measured by supplying it to positive bias(forward bias) and second: the cur- 
rent due to hole injection from Cu and Al was measured by reversing the polarity of the 
bias voltage, i.e., biasing Al and Cu electrodes positively(reverse bias). It is clear from Fig.l 
that J-V characteristics in case of ITO/CuPc/Cu display almost symmetric behavior in both 
cases(hole injection either from ITO, or from Cu electrodes). This is because in both cases 
there is small energy barriers of 0.05eV in case of ITO/MePc and 0.2eV in case of Cu/MePc 
interface, giving rise to space charge limited(SCL) bulk current when either ITO or Cu is 
positively biased. It has been shownj^, Ol, for Schottky energy barrier(SEB) less than about 
0.3-0.4eV, the current flow is due to SCL 3,33- But, in case of ITO/MePc/Al devices, J-V 
characteristics display diode like asymmetric behavior. Current density increases by almost 
five orders of magnitude by reversing the polarity of the bias. 

In several cases 0, Q, S|, the current transport mechanism has been explained in non- 
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crystalline organic semiconductor by extending the space charge limited current (SCLC) and 
trap charge limited current (TCLC) transport models in crystalline semiconductors! 
and J-V characteristics beyond Ohm's law was given by 



1G) and 



T/m+l 

Joc ^I (2) 
where m — 1 in case of SCLC and m = T /T in case of TCLC, T is the characteristic 
temperature of the traps and d is the thickness of the sample. We found in all our single 
layer devices with different electrodes(Al, Cu) and different MePcs(CuPc, ZnPc), J-V char- 
acteristics follow SCLC like behavior(J oc V 2 ) at low bias in forward bias condition, i.e. 
when ITO electrode is positively biased. It is appealing to attribute this behavior to SCLC, 
but as the bias increases the slope (m + 1) in log-log plot of J vs. V curves increases grad- 
ually from m = 1. This deviation can not be explained |i| as due to the presence of tr aps, 
because SCLC would not apply for low bias either. There are valid objections^, 3, (3 
against the use of transport theory(SCLC and TCLC) for crystalline systems in case of 
disordered systems, like OMS. Failure of the proposed SCLC and/or TCLC for crystalline 
11 1 in our devices and several other cases 0, can be argued for following 
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reasons, (i) As mentioned already, OMS does not have valence and conduction bands(as 
in case of crystalline semiconductor) due to lack of long range ordering in these materials. 
A distinguishing feature of molecular and polymeric organic semiconductors is that they 
are molecular glasses, composed of molecular entities held together by weak van der Walls 
interactions. They are characterized with HOMO and LUMO bands with GDOS. The tail 
states in HOMO/LUMO generally acts as traps and get filled with an application of bias due 
to the penetration of quasi- Fermi level inside the GDOS. (ii) SCLC and/or TCLC theories 
are applicable for crystalline systems and assume constant mobility. In case of molecular 
semiconductors, mobility is being affected by the disorders and the transport takes place 
by electrons hopping through MePc molecular sites having comparable energies within the 
GDOS resulting bulk limited behavior in the current conduction. In this case, mobility 
becomes a function of electric field due to the distribution in energy of hopping sites. Es- 
sentially, the electric field assists carriers in overcoming the potential barrier between sites 
with different energies. It has been found that the mobility fi of carriers in thin films of 
poly-phenylene vinvlenefPPV) P. lllH and tris, 8-hydroxyquinoline aluminum(Alq3)[14] over 
wide range of fields can be described according to universally observed PF relation(Eq.l). 



m 



iii) In several reports [7J, [a, |9|, thickness dependence of current has been used to corroborate 



the SCLC theory. But it has been shown |12| that current transport based on field dependent 
mobility due to hopping of charge carriers can account for similar thickness dependence. It 
is also possible to account for J oc V m relation without unphysical assumptions of TCLC 
theory. 

The temperature dependence of electrical conductivity associated with the hopping mo- 
tion in non-crystalline solids is explained by 



p(T) = p exp 



(3) 



'2o 
T 

where po is the temperature independent resistivity and T is the characteristic temperature, 



Qt — n-\- 1 

known as hopping exponent, D is the dimension, n is the exponent in the 
density of states g(E) in the vicinity of the Fermi energy Ep, which behaves like g(E) ~ 
(\E — Ep\) n . For an energy independent density of states(n = 0), the values a are 1/4 
and 1/3 in 3-D and 2-D, respectively. This is the Mott-Davis variable range hopping[15] 
mechanism. Shklovskii and Efrosjlfi| have shown that Coulomb interaction leads to a gap 
in g(E) around Ep and n = 1 in 2-D, i.e, g(E) ~ (\E — Ep\) and n = 2 in 3-D, i.e, 
g(E) ~ (\E — Ep\) 2 . Hence electron-electron interaction change the variable range hopping 
exponent to a — | in both 2-D and 3-D. Temperature dependence of resistivity p{T) over 
wide range of temperature for two thicknesses of CuPC is shown in Fig. 2. We have tried 
to fit the resistivity data with different a values. Fig.2 shows the ln(p) vs 1/T 1 ' 4 and also 
the ln(p) vs 1/T 1 / 2 plot. It is observed that the lOOnm structure shows better fitting with 
the 1/T 1 / 4 plot, whereas the 400nm structure shows better fitting with the 1/T 1 / 2 plot. 
For an improvement upon this eyeball approach, we adopt the method used by Hill ■ 
He introduces the quantity e a = d[ln(p)}/d(l/T) as the temperature dependent activation 
energy for variable range hopping conduction. Eq.3 can be expressed as 

In [^fjfy) = H*n) + (1 - a)lnT (4) 

One can determine a from e a vs T plot in log-log scale. We have plotted according to 
Eq.4, the straight line fit gives the value of a = 1/5 for 100nm(Fig.3) and a = 3/5 for 
400nm(Fig.3) CuPc layers. Since the value 1/5 is closer to 1/4 and 3/5 is closer to 1/2, the 
straight line fit in (j^j plots shows better fitting with a — | for lOOnm and with a = \ 
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for 400nm structures. Hence, it seems hopping mechanism may be in lower and higher 
thickness of CuPc layers. Moreover, hop ping mechanism may be altogether different from 
either Mott-Davis|lq or Efros-Shklovskii[l6j in OMS and the pre-exponential factor(p ) 
may depend on temperature. Fig. 4 shows the temperature dependent resistivity at various 
electric fields for ITO/CuPc(100nm)/Al structure. The linear dependence in lnp(T) vs 
1/T 1 / 4 upto 30K is observed. The electric field dependence on resistivity is studied for the 
Po at different electric fields. The straight line fit of ln(po) vs F 1 ! 2 plot, shown as inset 
in Fig. 4, signifies the universally observed Poole-Frenkle relation(/n(p) oc F 1 ^ 2 ) concerning 
charge carrier transport in OMS. 

To fit the experimental results, we have modified SCLC transport theory by simultane- 
ously solving the current density equation and Poission's equation with electric field depen- 
dent mobility (PF relation i.e. Eq.l) 



J(x) = p(x)ep [F(x), T] F(x) (5) 

e dF(x) , , , . 

where p(x) is the hole density at position x and p(F,T) is the hole mobility in MePc. 
Taking the boundary conditions (i) ohmic contacts at ITO/MePc interface(at x = 0) i.e. 
electric field vanishes(F = 0) at x — 0, and (ii) hole density is equivalent with the density 
of states(iV„) in HOMO of MePc, i.e, p(0) = N v = xl0 19 cm~ 3 , at x = 0. The excellent 
agreement between simulated J(V) results and experimental J-V characteristics for three 
different single layer deices based on CuPc and ZnPc over wide range of temperature is 
shown in Fig. 5. Following are the observations, (i) the simulated result follows J oc V m 
dependence, (ii) m is 2 at low electric fields and increases gradually from m = 2 at high 
electric fields and (iii) the value of m increases with temperature. Hence, the experimental 
data presented in Fig.2 , Fig. 3, Fig. 4 and Fig.5 alongwith simulation results establish the 
carrier transport process, which is thermally activated hopping with electric field dependent 
carrier mobility. 

There are three most important models for the explanation of particular temperature 
dependence of zero field mobility(p(0, T)) and the field activation energy(7(T)), appearing in 
universally observed PF relation(Eq.4.1) about charge carrier mobility. Those are, (a) Gill's 
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phenomenological model [J (b) Uncorrelated Gaussian disorder model [lj and (c) Correlated 
Gaussian disorder model 3]. Gill attempted to describe^] the experimental data with the 
following empirical form of charge carrier mobility 



fi(F, T) = no exp 
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where fio is the temperature independent mobility, A is the zero-field activation energy, B is 
a parameter, and T is the characteristic temperature of the material. In this case, /i(0,T) 
and 7(0, T) have 1/T dependence. Comparing this with the PF relation(Eq.l), fi(0,T) and 
7(0, T) can be expressed as 



MO, T) = fi exp 



k B Tj ,n ) k B \T T 



(8) 



The /j,(0,T) and j(T) at different temperatures, determined from theoretical fit to exper- 
imental J-V characteristics(shown in Fig. 5) are plotted in Fig. 6 according to Eq.8. From 
ln[(fj,(Q,T))] vs 1/T plot, the values of temperature independent mobility /io and the acti- 
vation energy A are calculated. The value of the constant B and the characteristics tem- 
perature To are calculated from the jCT) vs 1/T plot and given in Table-I. 

Bassler and co-workers proposedQ the uncorrelated Gaussian disorder model(UGDM) 
and provided support to PF behavior of mobility using Monte Carlo simulation. The UGDM 
describes the carrier transport as a biased random walk among the hopping sites with 
Gaussian-distributed random site energies. UGDM lead to following form of carrier mo- 
bility 



(jl(F, T) = fi exp 



r 



c 



a 



(9) 



where a = a/k B T, C & are constants, a is the Gaussian width of the hopping sites distri- 
bution. In this case, /x(0, T) and 7(0, T) have 1/T 2 dependence, instead of 1/T dependence 
as in case of Gill's empirical relation. Comparing this with the PF relation(Eq.l), MO, T) 
and 7(0, T) can be expressed as 



MO, T) = /i exp ( - - cr 2 ) , 7 (T) = C (a 2 



(10) 



Fig. 7 shows the temperature dependence of the /i(0,T) and 7(T) according to UGDM and 
the apparent linear dependence of ln/i(0) on 1/T 2 and 7(T) on 1/T 2 can be described by 



Eq.10. Values of temperature independent mobility //o, width of the energy spread ex, and 
the constants C and S are determined from the ln/i(0,T) vs. 1/T 2 and 7(T) vs. 1/T 2 
plots, and given in Table II. Although UGDM explains some features of experimental data 
and provides support for PF behavior of carrier mobility, several discrepancies emerge with 
uncorrelated description of Gaussian disorder model, which will be discussed later. The 
most important criticism against UGDM is its' inability to reproduce the PF behavior over 
wider range of electric field. Garstein and Conwellj^ first showed that a spatially correlated 
potential is required for the description of PF behavior of mobility for wider range of elec- 
tric field. Dunlop and co-workers [2j have shown that the interaction of charge carriers with 
permanent dipoles located on either dopant or host molecules give rise to PF behavior of mo- 
bility. Essentially, this correlated Gaussian disorder model(CGDM) is based on long-range 
correlation jjj between charge carriers and the molecular electric dipole, resulting random po- 
tential energy landscape with long-range spatial correlations, < U(0)U(r) >~ cr 2 a/r, where 
a is the minimal charge-dipole separation or the lattice constant of the molecular solid and 
a is the width of the Gaussian distribution and can be given by a = 2.35ep/ea 2 , where p is 
the dipole moment of the molecule. In this case, carrier mobility can be given by 



where A and T are parameters of the model. T characterizes the geometrical disorder. In this 
case, /x(0, T) has a 1/T 2 dependence (same as in UGDM) but 7(T) has a 1/T 3 / 2 dependence. 
Comparing this with the PF relation(Eq.l), /z(0, T) and 7(0, T) can be expressed as 



Fig. 8 shows the temperature dependence of /i(0, T) and j(T) according to CGDM, described 
by Eq.12. Values of the parameters [Iq, a, a and T are determined from the straight line 
plots of ln/i(0,T) vs. l/T 2 (shown in Fig.8) and 7(T) vs. l/T 3//2 (shown in Fig. 8) and given 
in Table III. 

Though the uncorrelated models(Gill's phenomenological relation and UCGM) can de- 
scribe the experimental data upto a some extent, but have certain limitations which are 
discussed below. There are several conceptual problems with Gill's phenomenological de- 
scription and the most important are (i) the obtained values of zero field and temperature 
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independent mobility (/xo) are overestimated by several order of magnitudes and don't match 
with the values obtained by independent measurements, (ii) the empirical Gill's formula 
lacks theoretical justification and the significance of different parameters(A, B and T ) are 
not well understood in the context of the physical properties of the system and (iii) the em- 
pirical formula for 7(T) predicts the negative field dependence for T > T , which deviates 
from the universal feature (PF) of carrier transport in these materials. 

Main feature of UGDM is the non-Arrehenious behavior of the temperature dependence 
of mobility, which is the consequence of hopping conduction in GDOS. But, there are several 
discrepancies in this description and those are (i) fitting is not as good as in case of CGDM, 

(ii) the intercept of j(T) vs. l/T 2 (shown in Fig. 7) gives positive value, but according to 
UGDM relation Eq.10, it should be negative, which is a fundamental problem with UGDM, 

(iii) it has been observed that there is no specific trend in the values of C, which is scaling 
parameter in the model and cannot be linked to a physical parameter of the system and 

(iv) its' inability to reproduce the PF behavior(which is a universal feature of charge carrier 
transport in these materials.) over wider range of electric field 

Here, we compare our experimental data with the uncorrelated(UGDM) and corre- 
lated(CGDM) theoretical models proposing the stretched exponential electric field depen- 
dence of /i. As already mentioned, at every temperature, /i(0,T) and j(T) are obtained 
from the experimental data presented in Fig.5, using the numerical workout and plotted in 
Fig. 8 which shows the temperature dependence of /x(0, T) and j(T) according to CGDM, 
described by Eq.12. Following observations have been made regarding the CGDM descrip- 
tion of the experimental data, (i) linear dependence of ln/i(0, T) on 1/T 2 and j(T) on 1/T 3 / 2 
is evident in Fig. 8 and straight line fit is excellent and better than that in UGDM, (ii) for all 
samples the intercept of j(T) vs. 1/T 3 / 2 gives negative values, according to CGDM(Eq.l2), 
(iii) the experimental value of T are 1.2 in CuPc and 1.4 in ZnPc, which are very close to the 
predicted value 1.97[? ], (iv) The width of the Gaussian distribution a is lOOmeV in CuPc 
and 120meV in ZnPc and similar values for a are commonly observed in other molecular 
solids, (v) The lattice constant (a) has been found to be 19 A in CuPc and 17A in ZnPc, 

n 

are in excellent agreement with the reported|21|] experimental values determined by X-Ray 
diffraction. 

The UCDM and CGDM share the common feature regarding the temperature depen- 
dence of //(0, T) due to similar distribution of hopping sites energies, but the temperature 
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dependence of field activation parameter 7 decides critically the importance of spatial corre- 
lation and decides the mechanism of carrier transport in organic molecular solids. Though 
CGDM explains the experimental data successfully, there are some limitations with this 
model. The positional disorder is not included in the present form of CGDM, but large 
values of disorder parameters a and T in NPPDA 3] , Alq3 22| and MePc signify the impor- 
tance of the positional disorder. In CGDM, it has been shown that long-range interaction 
between charge carriers and permanent dipole moments of doped molecules in polymers and 
host molecules leads to spatial correlation. However, it has been pointed out by Yu et. 
ai.j3] that the mechanism responsible for PF behavior in different conjugated polymers and 
molecules cannot be due to only charge-dipole interaction, because PF behavior has been 
universally observed in several doped and undoped polymers and molecules with or without 
permanent dipole moment. Hence, in addition to charge-dipole interaction there may be 
another mechanism responsible for spatial correlation, which is a fundamental requirement 
for PF behavior. Yu et. aI.[2J| have shown using first principle quantum chemical calcula- 
tion that the thermal fluctuations in the molecular geometry can lead to spatial correlation. 
It has been shown that the primary restoring force for the thermal fluctuation is steric and 
intermolecular, which lead to spatially correlated fluctuation in the energies of the localized 
states. Further experimental investigations are in progress to identify the exact origin of 
spatial correlation in molecules with and without dipole moments. 

In summary, we have studied the mechanism of carrier transport in OMS, which belongs 
to class of disordered molecular solids, by means of simple experiments. We have (i) shown 
that the carriers conducts through hopping; (ii) given evidences for the universal electric 
field dependence(/n(/i) oc VE) of carrier mobility and (iii) presented a comparison of the 
transport properties with uncorrelated and correlated disorder model. We have shown that 
E and T dependence of experimental data can be fitted and described within the correlated 
Gaussian disorder model, signifying the important role of both energetic and spatial disorder 
on the carrier transport in these materials. The origin of correlation has been discussed. We 
hope present analysis would lead to further experimental and theoretical investigations fo- 
cusing the exact origin and role of correlation on transport properties of disordered molecular 
solids. 
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Figure Captions 

Fig.l. Room temperature J-V characteristics for single layer hole only CuPc device with 
thickness of 200nm. Empty circles represent the data when ITO is positively bi- 
ased (Forward biased) and empty squares represent the data when ITO is negatively 
biased (reverse biased). The high rectification in (a) ITO/CuPc/Al and comparable 
currents in (b) ITO/CuPc/Cu ensures the condition to achieve the bulk limited current 
is when ITO is positively biased. 

Fig. 2 Temperature dependent resistivity of CuPc observed in ITO/CuPc/Al structure with 
lOOnm and 400nm of CuPc layers. Empty symbols are experimental data and solid 
lines are the straight line fit in the log(p) vs. 1/T 1 / 4 in (a) and vs 1/T 1 / 2 in (b). 

Fig. 3. Temperature dependent resistivity of CuPc in ITO/CuPc/Al structure plotted ac- 
cording to Eqn.4 for (a)100nm and (b)400nm CuPc layers. From the straight line 
fit (solid line) to the experimental data(symbols)), the values of a — 1/5 for lOOnm 
and a = 3/5 for 400nm structures were found out. 

Fig. 4. Temperature dependent conductivity of CuPc observed in ITO/CuPc(100nm)/Al 
structure. Empty circles are experimental data and solid lines are the straight line 
fit in the log(p) vs. 1/T 1 / 4 . Inset shows the electric field dependence of temperature 
independent resistivity (po) at different electric fields and the linear fit (solid line) shows 
the PF dependence of the carrier mobility. 

Fig. 5. J-V characteristics of bulk limited current in ITO/MePc/Al with ITO anode and Al 
cathode contacts at different temperatures starting at 320K and then at the interval 
of 30K. Experimental data are shown by the empty symbols and solid lines represent 
the simulated data, for (a) lOOnm CuPc layer (device A), (b) 200nm CuPc layer (device 
B), and (c) 200nm ZnPc layer (device C). 

Fig. 6. (a) Temperature dependence of zero field mobility /i(0, T) and (b) the field activation 
j(T) obtained from the theoretical fit to the experimental data shown in Fig. 5, are 
plotted according to the Gill's phenomenological observation (Eq. 8) for devices A, B, 
and C. The values of the parameters /x , A, B, and T are given in Table-I. 
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Fig. 7. Temperature dependence of zero field mobility /i(0, T) and the field activation 7(T) 
obtained from the theoretical fit to the experimental data shown in Fig. 5, are plotted 
according to the UGDM(Eq.lO) for devices A, B, and C. The values of the parameters 
fi , cr, C and E are given in Table-II. 

Fig. 8. Temperature dependence of zero field mobility /x(0, T) and the field activation 7(T) 
obtained from the theoretical fit to the experimental data shown in Fig. 5, are plotted 
according to the CGDM(Eq.l2) for devices A, B, and C. The values of the parameters 
/i , (cr), and a are given in Table-Ill. 

Table Captions 

Table. I Values of yU , A, B and T obtained using Gill's phenomenological relation for three 
different samples. 

Table. II Values of /x , cr, C and £ obtained using UGDM for three different samples. 
Table. Ill Values of /i , cr and a using CGDM for three different samples. 
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Table-I 



Device 




A 

A 


i > 

D 


J-o 




in — 

111 Vs 


in eV 


in eV(cm/V)2 


in K 


A 


2.76 x 10" 3 


0.44 


6.38 x 10" 4 


660 


B 


3.78 x 10- 3 


0.43 


5.94 x HT 4 


615 


C 


6.1 x 10" 4 


0.46 


5.28 x 10~ 4 


610 



Table-II 



Device 




a 


C 


£ 




(in cm 2 /V.sec) 


(in meV) 






A 


2.8 x lO" 8 


100 


6.4 x 10~ 5 


-3.1 


B 


3.0 x 10~ 8 


95 


6.9 x 10~ 5 


-3.18 


C 


1.7 x 10~ 8 


110 


4.5 x 10~ 5 


-4.24 



Table-III 



Device 




a 


a 




(in cm 2 /V.sec) 


(in meV) 


(in A) 


A 


2.8 x 10~ 8 


106 


20 


B 


3.0 x 10~ 8 


100 


19 


C 


1.7 x 10~ 8 


118 


17 
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10' 



CM 

E 


10" 9 


o 




■ 

Q. 




E 





< 




c 






10" 4 f 






o 






10- 5 [ 




10" 6 [ 




10" 7 = 




(b) 



g0 




o Forward 
□ Reverse 



2 4 
Bias(in Volt) 



Figure 1 




Figure 2 




Figure 3 




Figure 4 
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Figure 5 




2468 2468 

1000/T in K' 1 




Figure-7 




Figure 8 



